Nonrelativistic Chern-Simons Vortices 
on the Torus 



N. Akerblom 1 '* G. Cornelissen 2 '} 
G. Stavenga 3 '! and J.-W. van Holten 1 ^ 

1 Nikhef Theory Group, Amsterdam, The Netherlands 
2 Department of Mathematics, Utrecht University, The Netherlands 
3 Theoretical Physics Department, Fermilab, Batavia, IL, USA 

January 24, 2010 



Abstract 

A classification of all periodic self-dual static vortex solutions of the 
Jackiw-Pi model is given. Physically acceptable solutions of the Liouville 
equation are related to a class of functions which we term O-quasi-elliptic. 
This class includes, in particular, the elliptic functions and also contains a 
function previously investigated by Olesen. Some examples of solutions are 
studied numerically and we point out a peculiar phenomenon of lost vortex 
charge in the limit where the period lengths tend to infinity, that is, in the 
planar limit. 
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1 Introduction 

In this paper we study periodic, static vortex solutions of the Jackiw-Pi model 
[U [2]|^] This is a 2 + 1-dimensional nonrelativistic conformal field theory whose 
field content consists of a complex scalar field \& with non-linear-Schrodinger type 
action minimally coupled to a U(l) Chern-Simons gauge field A^. Let us begin by 
reviewing the elements of this model. 

1.1 The Jackiw-Pi model 

We take as our starting point the action [B] 

S[V, A„] = J dx° J d 2 x I - ^eij {AodiAj + AidjA + Ajd Ai) 

(1) 

1 a 2 ~\ 

+ ity*D ^ - - (D^)* • (D^) - — 1^| 4 |. 

1 For reviews see [3] 0] [5] . 
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Here E\ 2 = —£21 = +1, whilst = — ieA^ is the gauge covariant derivative 
and bold type indicates its spatial 2-vector part. We use the generic notation x° 
for the time coordinate, and apply the summation convention for repeated indices. 
In the following we define 



B = d 1 A 2 -d 2 A 1 , E i = d A i -d i A . 
The field equations derived from the action ([I]) then read 

B = ep, Ei = eeijJj, 

iD^ = --D 2 ^ + # 2 |^| 2 ^. 
2 



(2) 



(3) 



The chiral derivatives 

D ± = -L (A ± iD 2 ) , (4) 



satisfy the identities 



\ D 2 = D-D+ - "-B = D + D_ + e -B. (5) 



Using equations ([3]), the Schrodinger equation for ^ can then be written as 



iD ty = —D-D+'fy + yg 1 + —j |^| 2 ^ = -D + D^ + \^g 2 - —j |^| 2 *. (6) 

By a similar argument, the hamiltonian takes the form 

H = jd 2 x Q|Dvl/| 2 + y l^l 4 ) = / d 2 x (W| 2 + \(g 2 ± e 2 )\^ . (7) 

Hence there are two possibilities for constructing stationary zero-energy solutions: 

(J) = and g 2 + e 2 = 0, 

(8) 

(II) £Ltf = and g 2 - e 2 = 0. 

By stationary we mean that physical observables such as the particle density and 
current are time- independent. This is achieved by separating space and time vari- 
ables as 

* = e^, (9) 
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with p non-negative and time independent: d$p = 0. Any time dependence there- 
fore resides in the gauge-dependent phase oj. 

Substitution of either of the Ansatze (I) or (II) for \1/ and the coupling con- 
stants (e, g) simplifies the Schrodinger equation ^ to 

2 2 

iD^ = (eA - d to)^ = |^| 2 ^ => eA = d u T j P- (10) 

In addition, the real and imaginary parts of either condition D±ty = lead to the 
real equations 

eAi = diU ± Eijdj In y/p, (11) 

and as a result 

eB = eeijdiAj = =f A In y/p, A = d 2 + d 2 . (12) 
It follows directly that p satisfies one of the Liouville equations 
(/) D + ^ = Q => A In + e 2 p = 0, 

(II) D^ = =4> Aln^p- e 2 p = 0. 
The solutions of these equations are respectively of the form [7] 

4 \ff 



(13) 



(I) Pf 

(II) Pf 



e 2 (i + \f\ 2 r 

4 l/f 



(14) 



,2 (1 - i/| 2 ) 2 ' 

where /(z) is an analytic function of the complex coordinate 

z = x + iy, (15) 

and for physical reasons we make the hypothesis that / have at most isolated 
singularities (which then automatically are poles; see the discussion in Section 



13- 

Furthermore, boundedness of p requires \f\ 2 < 1 for case (II). This immedi- 
ately implies that there are no relevant non-trivial solutions in case (II). 

However, case (I) leads to a rich spectrum of vortex-type solutions, depending 
on the boundary conditions. For instance, if we take two-dimensional space to 
be the plane M 2 <H- C, it is necessary to require that at infinity p tends to zero 
sufficiently fast. The problem of writing down all static vortex solutions in this 
planar case was solved in a beautiful paper by Horvathy and Yera [8]. 
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For physical applications, e.g. in condensed matter systems, it is also of in- 
terest to study static vortex solutions in a finite volume with periodic boundary 
conditions; in that case one requires 

p(z + Ui) = p{z) (16) 

for given IR-linearly independent complex numbers U\, U2- This corresponds to 
studying the Jackiw-Pi model in the case where space is a two-dimensional flat 
torus C/(Zcji + ZfCo>2). Apparently, the first one to investigate this situation was 
Olesen, who gave a remarkable solution on a square torus [9]. 

If one thinks about how periodic boundary conditions are customarily employed 
in physics, the step from the plane to the torus seems innocuous. However, here 
this is not at all so. The change in topology, in fact, has a rather dramatic effect 
on the allowed rationalized vortex charg^] 

| vortex charge) (total magnetic flux) x e e 2 
q '~ 2vr x e ~ 2tt ~~ 2tt 

First of all, it can be shown that q, irrespective of whether we are working on the 
plane or the torus, must be a non-negative integer (see Appendix |B|. But, while 
from the classification of Horvathy and Yera it is immediate that on the plane the 
charge is always even: 

q = 2n, neN (vortex charge on the plane), (18) 

this is not the case for the Jackiw-Pi model on the torus. Indeed, Olesen's solution 
is an example of a static vortex on the torus with charge 

q = 1 (Olesen's solution on the torus). (19) 

Moreover, the correspondence between the solutions on the torus and those on the 
plane is somewhat involved. Olesen's solution, for instance, vanishes in the limit 
where the period lengths tend to infinity, and in Section [3] we give an example of 
a solution for which the charge q is halved as we pass from the torus to the plane! 
The adagium that the limit of a periodic solution, as the periods tend to infinity, 
gives a planar solution, fails dramatically in the case of Olesen's solution. 




1.2 Classification of vortex solutions 



In complex coordinates ( 15 ) the non-linear wave equation ( 13 ) for positive chirality 
fields of type (J) reads 

<9<91np + e 2 p = 0, (20) 



2 For a proof of equations (17) see [3J. 
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where d := {d\ — id 2 )/y/2 and d := (d\ + id 2 )/V^- 



The general solution (14, J) of this equation was discovered a long time ago 



by Liouville [7], who was led to the study of equations (14, ISzII) in connection 
with his researches on the theory of surfaces with constant intrinsic curvature^] (see 
also [HI [12]; in [13] solutions with vanishing boundary conditions on a rectangle 



were investigated). We shall frequently call pf defined by eq. (14, I) "the density 
associated with /." 

For our purposes, on physical grounds, we make the hypothesis that / is to have 
at most isolated singularities. This is because we want to interpret pf as a soliton 
(a vortex) and this interpretation is upset when / has a non-isolated singularity]^] 

We also demand that p be bounded. In fact, we impose the stronger condition 
that the total particle number in the spatial domain 

pd 2 x, (21) 

proportional to the total magnetic flux carried by vortices, is finite]^ In the case of 
a periodic p, boundedness automatically follows from continuity, as we can interpret 
p as living on a compact space (the torus), and in this case, boundedness is all we 



need for the integral (21) to make sense. For vortices on the plane, one obviously 
needs to supplant this with a suitable decay condition at infinity, see [8]. 
It can be shown: 

Lemma 1 (Horvathy-Yera [8]). Let pf be the density associated with a complex 
function f having at most isolated singularities. If pj is bounded, then the only 
possible singularities of f are poles, i.e. f is meromorphic in the plane. 

In the plane case this extends to infinity, so that / is a meromorphic function 
on the sphere, that is, a rational function: 

Theorem (Horvathy-Yera [8]). Let the density pj associated with f be a vortex 
solution of the Liouville equation on the plane. Then f is a rational function, i.e. 
there are polynomials P{z) and Q{z), such that 

f(z) = m 

f[) Q(zY 

3 On a surface of constant curvature the conformal factor of the metric in isothermal coordi- 
nates satisfies the Liouville equation; that is, if the metric is ds 2 — p(dx 2 + dy 2 ) with p > 0, 
then p satisfies equation (1 1 31) and e 2 is equal to the Gaussian curvature K of the surface. In this 



situation, the case K < is, of course, not excluded and corresponds to solution (14 77). It 



is known that equation (13 I) has no nowhere vanishing solution on the torus |10j . Thus, by 
necessity, all our torus solutions given below have zeros. 

4 Indeed, we may conjecture that if / has a non-isolated singularity, then its associated density 
Pf is unbounded. 

5 In the plane case the integral extends over R 2 , whereas in the periodic case it is taken over 
some elementary cell; say, the closure of the fundamental region: {t\UJi + t2^i \ < t\,ti < 1}. 
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Moreover, the converse is also true. 



In the case of the torus, Lemma [T] still holds (since it is a local statement), but 
boundedness of pf is automatic, and there is no corresponding statement about 
the behavior of / "at infinity." 

We now state the analogous classification in the case where p is periodic, or, 
as one could also say, lives on a torus: 



Theorem 1. Let p be a smooth periodic solution of the Liouville equation (|13j) with 
periods U\ and u 2 . It follows that p = Pf for some complex function f (Liouville, 
|7|/) meromorphic in the plane (Lemma^ which falls into one of the following two 
cases: 



Case A There are complex numbers p\, p 2 with \pi\ = 1, such that 

f(z + Wi) = Pif(z), 



(22) 



that is, f is an elliptic function of the second kind with multipliers pi of unit 
modulus. For the reader's convenience, we repeat the results of [Lflj for such 
functions in Appendix [C| 

Case B There are complex parameters z\,...,z n in the fundamental region of the 
lattice ZjUi + 7*uj 2 , and complex constants a$, . . . , a n , such that 



<p(z) 



1 



cp(z) 



where 



<p{z) 



k=l 



d k C, 



O(z) 



a(z - z ) 



(23) 



,C("i/2)« 



(24) 



with z = ^ + J J2k=i z k, and 



0{z) 



(25) 



for a suitable choice of parameters b, c, d, given in equations (65) and (66) 



Moreover, the converse is also true: If f falls into one of the two cases above, its 
associated density pf is a periodic solution of the Liouville equation. 

This result is derived in the next section. 

Remark on special functions. Our conventions for the appearing special functions 
are as follows: 
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• fpoji,w 2 indicates the Weierstrass p-function associated with the lattice = 

• C — Cwi,w 2 an d ° = a u>i,u2 are the Weierstrass zeta- and sigma-functions. 

The properties of these functions are given in many textbooks; see, for example, 
[15]. A word of caution: In the older literature, e.g. in the standard reference |16j . 
P = fpwi,aj 2 often denotes the Weierstrass p-function with half '-periods oji, uj 2 (and 
similarly for ( and a). 



2 Periodic vortices 

We now proceed to classify all periodic vortices on a given flat torus (Theorem [I]). 
To this end, let a lattice Q C C be given and suppose it is spanned by ui,u 2 , that 
is, 

0, = ZuJt + Ztu 2 . (26) 

As follows from our earlier discussion in Section the task is to find all smooth 
solutions p of the Liouville equation (13) such that 

p(z + oo) = p(z) for all u G fl (27) 

Suppose we are given such a p. Then, from [7] and Lemma [I] we know that there 
is a complex function /, meromorphic on the plane, such that 

4 |/'| 2 

p = Pf = ^(i+\jvr (28) 

where the prime ' denotes the derivative with respect to the complex variable z. 
Let w e (] be arbitrary and define the function 

g(z) :=f{z + u>). (29) 



From equation (27) and the fact that g'(z) = f'(z + u), it follows that 

p f (z) = p g {z) for all zeC. (30) 

In Appendix [A] we prove: 

Lemma 2. Let fi and f 2 be non-constant meromorphic functions on the plane 
and suppose that their associated densities pj 1 and pf 2 are equal: pf x = pf 2 . 



6 Another proof has been given by de Kok [T7] . 



8 



Then there exists a matrix 



such that 



7 



f 1 (z)= 1 -f 2 (z) 



a b 
c d 



G SU(2), 



a b 
c d 



Hz) :- 



af2{z) + b 



(31) 



a b 
c d 



cf 2 (z) + d 

Also, the converse is true 17^/ . even under the weaker hypothesis that 
U(2); that is, if f\ = V ■ f 2 for some V G U(2), then p h = p /2 Q 

Now, from Lemma|2]it follows that for any tu G fi, there is a matrix 7^ G SU(2), 
such that 

/(z + w) = </(*)= 7* •/(*)■ (32) 

This matrix is not unique in SU(2), but it is unique in PSU(2, C). We shall call a 
meromorphic function on the plane fi-quasi-elliptic if it satisfies condition (32). 

A trivial corollary to Lemma [2] is that p = pf is periodic with respect to the 
lattice Q = 7Lx)\ + if and only if there exists matrices 8 Ui G U(2), such that 
f(z + Ui) = 8 m ■ f(z) for 2 = 1,2. 

a b 

G SU(2) there is naturally associated a certain 



With every matrix 7 



c d 



transformation Ti^f) G PSU(2, C) from the Riemann sphere C to itself, namely 

az + b 



T( 7 ) : C -> C, 



z 1 — y 



cz + d 



Since, obviously, 



i) = ^(7^)^(7^) = T(^)T(^) for all u, u G Q, 



(33) 



(34) 



equation (32) tells us that any f2-quasi-elliptic function effects a group homomor- 
phism 

T-.n^G, u^T(>y u ), (35) 

from the lattice Q to some abelian subgroup G of PSU(2, C). We recall that 
PSU(2, C) = SO(3), the group of orientation preserving isometries of the sphere]^] 
Because Q = Zo>i + Zw 2 is a free module with generators U\ and u 2 , G is an 
abelian group with at most two generators T(j ui ) and T(j U2 ). 



r For M 



a 
7 S 



€ GL(2,C) and any complex function / we define M ■ f(z) 



af(z)+P 
~ff(z)+8 ■ 



For completeness we mention the elementary rule IX7172) = 2771)^(72) for all 71,72 G U(2). 
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Now, the important thing is that the converse part of Lemma [2] guarantees 
that any fi-quasi-elliptic function will also yield a periodic vortex solution of the 
Liouville equation. Therefore, the problem of finding all periodic vortex solutions is 
equivalent to writing down all f2-quasi-elliptic functions and this is directly related 
to classifying all abelian subgroups of PSU(2, C) with two generators. 

There are various ways to classify such subgroups. We will work in PSU(2, C) 
directly, and lift the two generators of the group from PSU(2,C) to SU(2). One 
can also use the isomorphism with SO (3), or consider rotations as quaternions. 
We shall comment on this later. 



By an earlier remark (immediately below equation (32)), we have the implica- 
tion 

T( 7 ) = T(7) 7 = ±7 for all 7, 7 G SU(2). (36) 
Then, since the generators T( 7cJ1 ) and T{^ W2 ) of G commute, 

T(7 Wl )T( 7w2 ) = T( 7t , 2 )T(7 Wl ), 

it is easy to see that 7cJl and j W2 either commute or anticommute. We will refer 
to these cases as Case A and Case B, respectively and treat them in turn in the 
following two sections. 



2.1 Case A: The matrices 7^ and 7 W2 commute 



Since 7cJl and 7 W2 commute, they can simultaneously be put into diagonal form. 
More precisely, there exists a matrix U G SU(2), such that 



1 1 Jim 



u 



1,2), 



(37) 



where the are complex numbers of unit modulus: 
Let / be fi-quasi-elliptic and define the function 

g(z) = U-f(z). 



1. 



(3* 



It follows that 

g(z + Ui) = iM9(z) (i = 1,2), (39) 

i.e. the function g is a so-called elliptic function of the second kind. There exists 
a complete classification of all such functions (cf. Appendix [C]). Thus, / will be of 
the form 

f = U^-g (40) 

with g some elliptic function of the second kind, and, by Lemma |2j the densities 
associated with these functions are the same: 



PS = Pa- 



(41) 
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Conversely, if g is a quasi-elliptic function of the second kind with multipliers 
Hi satisfying = 1, then its associated density p g is periodic. Indeed, for any 
such function g there are matrices 



Hi 



GSU(2) (i = l,2), 



with 



(42) 
(43) 



g(z + u)i) = 7 Wi • g{z) (i = l,2), 
and the claim immediately follows from the corollary to Lemma [2} 

2.2 Case B: The matrices 7^ and 7 W2 anticommute 

If our matrices j ui and 7 W2 anticommute, we can diagonalize one of them and 
anti-diagonalize the other. Specifically, there is a matrix U G SU(2), such that 



-% 
% 



7«* = ^ 

for some complex A with |A| = 1. Now put 

M : 

whence 



-A 

A- 1 



U. 



(44) 



1 

iX 



which is to say 



-i 
i 



7u* = ^ 



-i 
i 



V, 7^ = ^ 



z 

1 



Mt/, 



i 

1 



V 



(45) 



(46) 



(47) 



for some V = MU G U(2). 

Let us briefly digress to remark on the subgroup G of PSU(2, C) generated by 
T(7 UX ) and T(^ 2 ). 

If we define 



a := T(7 W1 ) : z i-» 



-i 
z 



and 



c := a o b : z >-)■ 



& := ^(7^2) : 2 !-)■ 



1 
-1 



i 

1 



z, (48) 



(49) 
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we get the composition table of the famous Vierergruppe V = Z 2 x Z 2 : 



o 


1 


a 


b 


c 


1 


1 


a 


b 


c 


a 


a 


1 


c 


b 


b 


b 


c 


1 


a 


c 


c 


b 


a 


1 



(50) 



where 1 denotes the identity transformation z H- z. Our subgroup G is isomorphic 

to Z 2 x Z 2 ! 

Coming back to our classification problem, it follows that any fi-quasi-elliptic 
function / is of the form 

f = V*-g, (51) 
where g is a function meromorphic in the plane satisfying 

g(z + ui) = -g(z), g(z + w 2 ) = ^/g(z). (52) 

Conversely, from the corollary to Lemma [2] it is plain that the density pj associated 
with any such / is periodic, for there are matrices Mi,M 2 e U(2), such that 
f(z + oji) = Mi ■ f(z) for i — 1,2. Moreover, pf = p g . 

We now proceed to classify all meromorphic functions in the plane which satisfy 
the period condition (52). Suppose go(z) is some such function satisfying equation 

(53) 



(52) and let g(z) be any other such function. Put 

f(z) := g(z)/g (z). 



Then 

If we define 
with 

it follows that 



f(z + u 1 ) = f(z), f(z + u 2 ) = l/f(z) 
<p(z) := IP ■ f{z) 

U :-- 



-1 1 
1 1 



(p(z + ui) = (p{z), ip(z + u 2 ) = -<p(z); 

therefore, <f(z) is some multiplicative quasi-elliptic function with p,\ 
From Appendix O we find that there are complex constants 



1, 



and parameters 



ciq, . . . , a n G C, 
Z n E {hUi + t 2 id2 I o < tl, t 2 < 1} 



(54) 
(55) 

(56) 
(57) 

: -1. 

(58) 
(59) 
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in the fundamental domain of the lattice f2 = Zc^i + Zw 2 , such that 



tp(z) 



k=l 



a(z-z ) n cC(bn/2)l 



where 



i n 

z o — ~ r - ) z k . 

In n ^-^ 



k=l 



Therefore, g is of the form 



g{z) = [U - <p(z)] 9o(z) = -^^\go(z) 



(60) 



(61) 



(62) 



Conversely, any such function g satisfies the conditions (52). 

It remains to give some g$ satisfying equation (52). Inspired by Olesen's special 
solution j!) . we make the Ansatz 



9o(z) = 0{z) 



Cp2u 1 ,2uj 2 { z ) + d' 



(63) 



We have the general formulas 

(ei - e 2 )(ei - e 3 ) 



p(z + ui) = ei + 



/ i \ . (e 2 -ei)(e 2 -e 3 ) 

p(z + w 2 )=e 2 + ^ , (64) 

- e 2 



with p = p2u 1 ,2u 2 , ei := p(u>i), e 2 := p(o> 2 ), and e 3 := -(ei + e 2 ). Using these 
formulas and demanding that g satisfy (52), we can choose the parameters b, c, 
and d in our Ansatz (63) appropriately. With the help of a computer algebra 

c(— 2ei + e 2 - c 2 e 2 ) 



system (Mathematica) we have found that 

-e 2 + c 2 (-2 ei + e 2 ) 
1 + c 2 

with 



1 + c 2 



-3ex + -e 2 )(2ei + e 2 ) 



(65) 



(66) 



ei + 2e 2 

will do, as long as ei + 2e 2 7^ 00 Indeed, e% + 2e 2 = only in the limit where our 
torus degenerates into a cylinder and this is excluded. This concludes our proof of 
Theorem [TJ 



9 It turns out to be immaterial which branches we choose for the square roots. In this sense, 
the choice of parameters is essentially unique. 
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2.3 The abstract underlying group 

We now explain how to refine our classification from a different perspective, using 
the isomorphism PSU(2, C) = SO(3) = H 1 with the different model groups of space 
rotations, and unit quaternions H 1 . Let us denote by G the subgroup (in any of 
these models) generated by j ui and 7 W2 . Then G is an abelian group of rotations, 
which is intrinsically attached to the vortex solutions of the torus Jackiw-Pi model. 
We call the abstract isomorphism type of this group the type of the vortex solution. 

We denote by Q the set of rational numbers. As usual, we call a real num- 
ber irrational if it is not rational. We call two real numbers linear dependent 
over Q (abbreviated "LD") if one is a rational multiple of the other (and linear 
independent otherwise). 

Suppose our rotations are around the same axis, one through an angle 27r#, the 
other through an angle 2ir8'. If one of 9 and 6', say 9, is rational with denominator 
m, then its associated rotation generates a cyclic subgroup Z m of G of order 
m. If then 9' is irrational, we find that G = Z m x Z (where it is possible that 
m = 1, in which case G is infinite cyclic: G = Z). If both 9 and 9' are rational 
with denominators m and n, say, then G is a cyclic group of order the least 
common multiple lcm(m,n) of m and n, that is, G = Zi cm ( mj „), a finite cyclic 
group (possibly trivial, which corresponds to genuinely elliptic functions). Finally, 
if 9 and 9' are both irrational and linearly independent over Q, the corresponding 
rotations generate a group G = Z x Z, but if they are linearly dependent over Q, 
they generate a group G = Z. 

Suppose now that G consists of two commuting rotations around different axes. 
It is easy to show (e.g., using the unit quaternion picture, in which a rotation 
around an axis v = (vi,V2,vz) through an angle 29 is represented by cos 9 + 
sin 9{v\i + V23 + v$k)) that the only pair of commuting rotations are two rotations 
of 180° around two orthogonal axes, and then, abstractly, the group G is the 
Vierergruppe. Also, up to an isometry of space, we can assume that the axes are 
in a fixed position, so this group G can be conjugated in SO (3) into standard form. 

Thus, we see that Case A corresponds to rotations around the same axis, 
whereas Case B corresponds to the Vierergruppe of two rotations around two 
different axes. 



We have summarized the preceding discussion in Table 2^3 In this table, we 
denote by ord(/i) the multiplicative order of a complex number \i in C*, i.e., the 
smallest positive integer N for which /i N = 1 (and we put ord(/x) = 00 if no such 
integer exists). We call two complex numbers \i\ and ^2 multiplicatively dependent 
(abreviated "MD") if there exist integers N± and N 2 such that fi^ 1 = fi^ 2 - We 
denote a space rotation around an axis v through an angle 9 by R${9). Note again 
that in this table m and n are integers, so Zi cm ( m>n ) can be the trivial group (if 
m — n — 1), and Z m x Z can be an infinite cyclic group = Z (if m — 1). 
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in SU(2) 


in SO(3) 




type 


Case A: commuting 


Same rotation axes 






• ord(/ix) = m and ord(/x 2 ) = n 


(Rv{2<Klm),Rz{2<Kln)) 




^lcm(m,n) 


• ord( / u 1 ) = m and ord(/z 2 ) = oo 


(Rz(27r/m),Rv(27r9)), 9 <£Q 






• ord(/ii) = ord(/i 2 ) = oo MD 


(Rt(2ir6),Rt(2ii6')), 9,9' £Q 


LD 


z 


• ord(/ii) = ord(// 2 ) = oo not MD 


(Rs{2Tr9),Rz{2ir9')), 9,9' ^ <Q 


not LD 


Z x Z 


Case B: anticommuting 


Orthogonal rotation axes 

(Rt(7r),Rti(<K)) (v±w) 




Z 2 x Z 2 



Table 1: Possible "types" of vortex solutions on the torus (m, n are integers). 

3 Examples 

3.1 Flux loss and flux conservation for elliptic function so- 
lutions 

A brief glance at Theorem [T] will convince the reader that, in particular, the den- 
sities associated with elliptic functions furnish examples of periodic vortices (take 
/^i, /X2 = 1 in Case A). The type of these solutions is trivial. 

A function / is elliptic with respect to the lattice Q = Zo>i + Zo; 2 precisely if 
it can be expressed as 

f(z) = R 1 (p(z)) + p'(z)R 2 (p(z)), (67) 

where Ri, i? 2 are rational functions and p = p^ xW - 

It is easy to see that if we put u>i —> tuji (i = 1, 2) and take the limit t —> +oo, 
then (compare (20], pp. 85 ff.) 

f(z)^R 1 (z- 2 )-2z~ 3 R 2 (z~ 2 ). (68) 

That is, in the limit where we remove the periodic boundary conditions (the planar 
limit), / tends to a rational function. Since any rational function can be written 



in the form (68) for some rational function R 2 , any rational function can arise in 
this way as the limit of an elliptic function. Thus, in this way we obtain all static 
vortices on the plane. 

An elliptic solution with flux loss. Let pf t (t > 0) be the density associated 
with the function 

/.(*) := ^44- (69) 

(We are dealing with the torus C/ (Zt+Zit).) Figure [l] shows a plot of this density 
for t — 1. Numerical integration suggests that for the rationalized charge oWus 
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Figure 1: Plot of the density (in units of 1/e 2 ) associated with the function 
Pt it( z ) / Pt,u( z ) for t = 1 in the cell 0.7 < x < 1.7, 0.3 < y < 1.3. Large val- 
ues of the density have been clipped. 



associated with this solution (t > finite) we have 



- £l 

Q'torus „ 

Z7T 



^ Pft d 2 x = 4 (F:= [0,t] x [0,t]). (70) 



Now, the planar limit of f t is 



/*(*)-> — fort^+oo, (71) 

z 

and it is well known that the charge associated with this is 

e 2 r 2 x 

Qplane = ^~ / Pz^-2/z d X = 2 = - ^torus- (72) 

2ti J R 2 2 

We therefore have the surprising result that, in passing from the torus to the plane, 
some charge of a vortex can get lost. 



An elliptic solution with no flux loss. That this need not always happen is 
shown by the example of the density associated with g(z) = p t ,u{z)- Here, the 
charge in the planar limit is the same as on the torus, namely = 4. 
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3.2 Relatives of Olesen's solution 



In [9] Olesen investigated a periodic vortex with charge q = 1. In our language, 
this solution is associated with the function O (z) (equation (63)) for the square 
lattice Q = 'Lt + 'Lit with t > 0. In Figure [2] we have plotted this density for t = 1. 




Figure 2: Plot of Olesen's density (in units of 1/e 2 ) on a 2 x 2 grid of cells where 
the fundamental domain is [0, 1) x [0, 1). 

We can also look at the density associated with O (z) on arbitrary tori C/ (Zcui+ 
Zoj 2 )- For instance, Figure [3] shows the density for a sequence of lattices Cl = 
Z + 1M, where successively t = .5, .75, 1. Note how the drempel-like structure 10 
deforms to the lump of Figure [2] as the rectangle approaches a square. From nu- 
merical integration we know that all these vortices have charge q = 1 and the same 
appears to be true for tori where the fundamental region is a true parallelogram. 

What is the planar limit of the density associated with O (z)? It is easy to 
see that for a square fundamental region, O (z) approaches a constant as the pe- 
riod lengths tend to positive infinity, that is, in this limit, the associated density 
approaches 0. It appears likely that the same is also true for more general funda- 
mental regions. 

io "Drempel" is a Dutch word which, amongst other things, denotes a speed bump. 
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4 Summary and discussion 



In this paper we have studied the Jackiw-Pi model with periodic boundary con- 
ditions, which amounts to solving the Liouville equation on the torus. Physically, 
these solutions describe a two-dimensional periodic lattice of charged vortices with 
quantized magnetic flux. As first discussed in [fl [2] the existence of vortex solu- 
tions requires a delicate tuning of the coupling parameters: the electric charge 
and the strength of the self-interaction. Surprisingly, it seems that this tuning is 
not destroyed by quantum fluctuations [2TJ El [18] ; on the contrary, the tuning is 
precisely the condition for which the /3-functions of the model vanish and there is 
no scale-dependence of the parameters, at least at one-loop order. 

On the torus, the spectrum of fluxes of the vortices differs from the planar 
case; it is richer in that it allows both odd and even integer fluxes. This is possible 
because periodic functions on the plane do not vanish at infinity, as required for 
the solutions on the infinite plane. However, it also implies that the limit of the 
torus to the infinite plane is singular and can change the flux associated with a 
certain solution. We have presented explicit examples of this phenomenon. This 
observation may be relevant also in other field theories with soliton solutions, e.g. 
the Skyrme model as an effective theory for the bound states in QCD. 

It is amusing to note that our physical classification of vortices on the torus 
has a purely mathematical consequence having to do with the geometrical content 
of the Liouville equation: We can interpret our density p as the conformal factor 
of a metric on a punctured torus, with punctures exactly at the zeros of p. Our 
classification theorem then gives all sufficiently smooth metrics of constant Gaus- 
sian curvature K = e 2 > on punctured tori in explicit formj^] From our physical 
arguments in Appendix [B] it also follows that the properly normalized integral (17) 
of the conformal factor over the torus is always a non-negative integer. 

In reference |22j a topological interpretation of the charge of vortex solutions 
on the plane was given. It would clearly be interesting to obtain an analogous 
interpretation for the theory on the torus and we believe that the remarks in 
Appendix [B] could constitute the first steps in that direction. 
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A Proof of Lemma E 



Here we prove Lemma [2] of Section [2] We only need to supply the proof of the 
"^"-direction; for the "^"-direction see [HI HE]- For clarity, let us repeat the 
statement (in slightly altered notation): 



Lemma 2 ("=>"). Let f and f be non-constant meromorphic functions on the 
plane and suppose that their associated densities pf and pj are equal: pf = pj, 
where 

4 \f{z)\ 2 

(l + l/c*)! 2 ) 3 ' 



Pf( z ) 



(73) 



and analogously for pj. 

Then there exists a matrix 



such that 



7 



f(z) = 7 • f(z) 



a b 
c d 



E SU(2), 



a b 
c d 



af(z) + b 
cf(z) + d- 



(74) 



Proof. Stereographic projection it : S 2 — > C w gives a bijection between the sphere 
S 2 and the extended complex w-plane C w . In this way, the round metric on the 
sphere, ds 2 s2 , induces a distance function du on C w , for which the distance between 
two points wi, W2 G C w is given by 



d u (w 1 ,w 2 ) 



inf 

r 



\r'(t)\ 



o i + |rw 



dt, 



(75) 



where the infimum is over all curves T : [0, 1] — > C w with T(0) = Wi, T(l) = w 2 . 
The orientation preserving isometry group of the sphere, SO (3), is mapped by 7r to 
the orientation preserving isometries of C w equipped with the distance du, which 
isPSU(2,C). 

For a meromorphic function on the plane C z , define a quasi-distance df by 



df(z!, z 2 ) := du(f(zi), f{z 2 )). 



(76) 



(We call this a quasi-distance since, although it is positive and satisfies the triangle 
inequality, it is degenerate in the sense that points Z\, z 2 at distance zero are not 
necessarily equal, but rather satisfy f(z\) = f{z 2 ).) 
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The hypothesis of the theorem concerning equality of densities implies that for 
every Zi, z 2 G C z , we have 



df(zi, z 2 ) = df(zi, z 2 ). (77) 
We now define a map i : C w — > C w by 

l(w) := /(/-»)■ (78) 
First of all, this is well-defined. Indeed, if f(zx) = f(z 2 ) ='■ w, then the definition 



(76) implies that df(zi, z 2 ) = 0. Further, equation ( |77| ) implies that df(zi, z 2 ) = 0, 
and, again by definition (76), we obtain f(zi) = f(z 2 ). Our claim is that i(w) is 
an isometry of C ra equipped with d\j. 

It is surjective, since / and / are not constant. Indeed, for any two points 

W\,w 2 G C w , we have 

d v (L(wx),L(wJ)) = dj{f~ 1 {w 1 )J~ 1 {w 2 )) = d f {f~ 1 {w 1 )J~ 1 {w 2 )) =du(w u w 2 ). 

(79) 

Also, l is orientation-preserving since / is meromorphic. 

Hence, / = T(f) for some orientation preserving isometry T of C w , that is 
T G PSU(2, C), whence there is a matrix 



7 

such that 



a b 
c d 



G SU(2), (80) 



/ = 7 ' /• (81) 

□ 

Remark. It is clear that this lemma can be used for determining the precise struc- 
ture of the moduli space of self-dual static vortices of the Jackiw-Pi model on 
the plane. For, according to Horvathy and Yera [8], any such vortex with flux 
$ = 4nN/e is given by a density pf, where / is a rational function 

f{z) = wr de z p<de zQ = N - ( 82 ) 

Therefore, every such solution has 4N moduli but, obviously, they are not all 
independent. Rather, by our result, the moduli space is some kind of quotient 

C 2A 7PSU(2,C). 

The invariant theory of PSU(2,C) is well-studied, see e.g. |23j. We leave the 
problem of working out the physical implications in detail for the future [21]. 
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B Quantization of flux 



We comment here on the quantization of flux of static vortex solutions of the 
Jackiw-Pi model. 

For the theory on the plane, this quantization is best seen a posteriori from 
the results of Horvathy and Yera [8J. For the time being, an analogous result on 
the torus is, however, not available [24]. That is, given a solution from the classi- 
fication Theorem [T] we cannot say at the moment, without resorting to numerical 
integration, what its associated flux is. 

Therefore, we now proceed to give a more general argument supporting the 
claim that the flux is also quantized in the torus case. 

The boundary conditions of the Jackiw-Pi model on a spacetime of the form 
R x T 2 , where T 2 = C/Q for some lattice Q, are somewhat subtle. Naively, one 
would write the gauge potential A as a 1-form on the torus, which would lead to 



B = dA= A = 0, (83) 

T 2 JT 2 JdT 2 =® 

in contradiction to the solutions with a non-vanishing magnetic flux. The resolu- 
tion to this puzzle is of course analogous to the Dirac monopole, where we need 
multiple gauge patches to describe the solution; in other words, A in reality is a 
section of a bundle. 

However, because we are dealing with a torus, we can also pull back the gauge 
connection to the plane, where the gauge potential can be written as a 1-form. 
The boundary conditions are then implemented by periodicity of the fields p, E, 
and B, which translates to the equations 

A{x + u l ) = A(x) + d6i{x), ' 

where Q = 7Lu\ + Zw2, that is, our lattice is generated by U\ and U2- 

Now we can use gauge transformations in the plane to set the phase Qi to 
zero, and then we are left with a single phase Q\. It is easy to show that under 
translation by U2 we have 

e i$l(x) _ e i6±(x+u 2 ) ^gg^ 

and thus 6\{x + U2) = 6\(x) + 27m. This means that the total magnetic flux 
through the torus is 



B= dA= A, (86) 

F JF JdF 

where 

F:={t 1 u} 1 + hv 2 \0<t 1 ,t 2 <l} (87) 
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is the fundamental domain of the torus in the plane. The boundary integral is the 
integral along the parallelogram where the two sides in the direction of ooi cancel, 
due to periodicity of A in w 2 - However, the sides in the direction of u 2 do not 
cancel, due to the non-periodicity caused by Q\. The difference between the two 
sides is given by 

r r^2 

/ A— d6 1 = 2nn. (88) 

J OF JO 

Therefore, the total magnetic flux is quantized in units of 2rr. The topology of the 
principal U(l) gauge bundle over the torus is that of a twisted 3-torus with twist 
n. 



C Elliptic functions of the second kind 

For easy reference we repeat here the results of [14J, p. 154 concerning elliptic 
functions of the second kind (=multiplicative quasi-elliptic functions) specialized 
to the needs of the present paper (see also [25j 126]). 

Definition. Let Q = Zcji + Zcu2 be a lattice. A function / which is meromorphic 
in the plane is said to be an elliptic function of the second kind with multipliers of 
unit modulus, if there exist complex numbers /ii,//2, with I//2I = 1, such that 



/(* + £*)= (i = l,2). 



59) 



Theorem (Lu [E]). A function f which is meromorphic in the plane is an elliptic 
function of the second kind with multipliers \i\,\i 2 of unit modulus if and only if 
there are complex constants 

a ,...,a n eC, (90) 

and parameters 

zt,...,z n e {ti^i + t 2 u 2 I < t u t 2 < 1}, (91) 



such that 



where 



f{z) 



n d k ( 



k=l 



A 



' dz k 
1 

7TZ 



aiz - z ) 



(12m -71 m), 



and 



z 



1 1 n 

2^i^-^ ) + n^ 



Zk- 



(92) 

(93) 
(94) 



k=l 



Here, ^ := Cn^aM^) and 7* 
chosen arbitrarily.) 



log/Xj (i = 1,2). (The branch o/log/i, can be 
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(a) t = .5 



(b) t = .75 




(c) t = 1 



Figure 3: Density (in units of 1/e 2 ) associated with the function O(z) on a sequence 
of rectangular tori with lattice = Z + lAt for: [(a) t = .5, (b) t = .75, and [(c 
i = 1. 
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